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We study excited states of the nucleon in quenched lattice QCD with the spectral analysis using the maximum 
entropy method. Our simulations are performed on three lattice sizes 16^ X 32, 24^ X 32 and 32^ x 32, at = 6.0 
to address the finite volume issue. We find a significant finite volume effect on the mass of the Roper resonance 
for light quark masses. After removing this systematic error, its mass becomes considerably reduced toward the 
direction to solve the level order puzzle between the Roper resonance A^'(1440) and the negative-parity nucleon 
7V*(1535). 



A particular interest in the excited nucleon 
spectra is the level order of the positive-parity 
excited nucleon A^'(1440) (the Roper resonance) 
and the negative-parity nucleon 7V*(1535). This 
pattern of the level order is also seen universally 
in the excited states of flavor octet and decuplet 
baryons. The quark confining models such as the 
harmonic oscillator quark model or the MIT bag 
model have difficulties in reproducing the Roper 
as the first excited state of the nucleon ^ . 

First systematic lattice QCD calculation for 
the nucleon excited states in both parity chan- 
nels shows that the wrong ordering between N' 
and N* also takes place at least for heavy-quark 
masses After this work, many lattice calcu- 
lations confirmed this puzzle, which are summa- 
rized in Ref. 

Recently, it was pointed out that the finite lat- 
tice volume affects the structure of the Roper res- 
onance considerably for light quark masses |3I4) . 
If the "wave function" of a quark inside baryon 
is squeezed due to the small volume, the kinetic 
energy of internal quarks increases and the total 
energy of the bound state should be pushed up. 
Such effect is expected to become serious for ra- 
dially excited states as compared to the ground 
state. In the present report, we show that the 
Roper resonance particularly receives large finite 
volume effect: its mass after removing the finite 
volume effect moves toward the direction to solve 



the level order puzzle. 

In the numerical simulations, we generate 
quenched QCD configurations with standard 
single-plaquette Wilson action mainly at /? = 6.0 
(fl-i ~ 2GeV) for three diflFcrcnt lattices, 16^ x 32 
(444 conf.), 243 x 32 (350 conf.), and 32^ x 32 
(200 conf.) to study finite volume effect. The 
quark propagators are computed using the Wil- 
son fermion action at four values of the hop- 
ping parameter, which cover the range M^/Mp ~ 
0.69—0.92. We adopt the conventional interpolat- 
ing operator eabc{'u^Cj5db)uc to create the posi- 
tive parity states of the nucleon. 

To extract the mass of the excited nucleon, we 
utilize the maximum entropy method (MEM) 
which can reconstruct the spectral function (SPF) 
A{uj) from given Monte Carlo data of two point 
hadron correlator G{t) — J dijjA{uj) ex-p{—ujt). 
Details of the MEM analysis can be found in |4I5| . 

First we show dimensionless SPF of the nucleon 
for the largest lattice size at the lightest quark 
mass as a typical example in Fig^ There are two 
sharp peaks and two large bumps. The crosses 
on each peak or bump represent the statistical 
significance of SPF obtained by the MEM 0. 
We remark that the first two peaks are phys- 
ical states but two large bumps are unphysical 
states. This is confirmed by the observation that 
those bumps appear at the same dimensionless 
frequency Lua in simulations with different lattice 
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Figure 1. The dimensionless SPF ^(0-1) — A{lu)uj^ 
in the nucleon channel as function of the fre- 
quency Lu in lattice units. 

spacings. It indicates that those states become in- 
finitely heavy and decouple from physical states 
in the continuum limit. The similar observation 
in the mesonic case has been reported in Ref. |^. 

In Fig.[21we plot the masses of the ground state 
(open circles) and the first excited state (filled 
circles) , which correspond to the peak positions of 
first two peaks in SPF, as functions of the spatial 
lattice size L. The quoted errors are estimated 
by the jackknife method. Observed finite volume 
effect is clearly sizable for light quarks (the right 
panel) and its effect is larger for the excited state 
than the ground state. Dashed horizontal lines 
in Fig[21 express the masses in the infinite volume 
limit and their errors guided by a formula = 
Moo+cL~^ j7|- For the lightest quark mass in our 
calculation, L should be as large as 30 to remove 
the finite volume effect. This corresponds to La w 
3.0 fm in the physical unit. 

To see whether the level switching between N' 
(the Roper) and N* (the negative parity nucleon) 
may possibly occur, the squared masses of these 
resonances are plotted in FigOl as a function of 
the pion mass squared. The filled circles arc our 
results of the Roper in the infinite volume limit, 
while the open squares and triangles are the data 
at the same lattice spacing for the negative par- 
ity nucleon taken from Refs. j2l8j for La Ki 1.5 
fm and 2.2 fm. Although the negative parity 
state is lighter for heavy quark masses, there is a 
tendency that the N'-N* level splitting becomes 
smaller or even inverted as the quark mass de- 
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Figure 2. Finite volume effect on the positive- 
parity nucleon channel. Masses of the ground 
state (open symbols) and the first excited state 
(filled symbols) for three lightest quark masses as 
a function of the spatial lattice size L. 

creases. Note here that the finite volume effect is 
not significant for the negative parity state jH] as 
can be seen by the comparison among the open 
symbols in FiglHl This is in contrast to the signif- 
icant finite volume effect on the Roper resonance 
shown in Fig[21and in Refs. |3I4| . 

Taking a simple linear extrapolation toward the 
chiral limit in FigEl we find A/jv = 0.73(14) in 
lattice units. Our N' mass is evidently smaller 
than the previously published results for the N* 
at the same lattice spacing; Mjv* = 0.85(5) [2] 
from domain wall fermions and AI^* — 0.89(2) [S] 
from clover fermions. 

Next we demonstrate the importance of the 
large spatial volume in a different way. We cal- 
culated the quark mass dependence of the mass 
splitting between the ground state and radially 
excited state. From the phenomenological point 
of view, such a mass splitting is almost indepen- 
dent of the quark mass as pointed out in Ref. P] . 
Shown in Figgis our lattice results of the N-N' 
mass splitting for the largest size lattice La = 3 
fm. As is evident from the figure, the mass split- 
tings are independent of the quark mass within 
the statistical error bars. This result is different 
from those shown in Ref. [HJ. 

In summary, we studied the excited states of 



3 



o (1/2) from DWF, La~1.5fm 
V (1/2)' from Clover, La~1 .5frr 
a (1/2)' from Clover, La~2.2fm 



(1/2) 1 St excited state (this work) 
(1/2)"*" ground state (this work) 
J I l_ 



_L_ 



0.0 



0.1 



0.2 



0.3 



0.4 



0.5 



(aM;i)' 




(aM;:)'^ 

Figure 4. The splitting of 1st and 2nd peaks 
of the spectral function as a function of the pion 
mass squared. The corresponding experimental 
values for M^^^^^^^ - M^^^^^j and for the strange 
sector, M^(^^^„, -^53(1190) and A/^^^^^^, - M^^^^^^, , 
are marked as stars. 



Figure 3. Ground and excited nuclcon spectra 
as a function of the the pion mass squared. The 
corresponding experimental values for iV(940), 
7V'(1440) and iV*(1535) are marked with lower, 
middle and upper stars. 

the nucleon with the use of MEM. On the basis of 
a systematic analysis with three different lattice 
volumes {La = 1.5, 2.2 and 3.0 fm), we confirmed 
our previous finding that the Roper resonance re- 
ceives large finite volume effect for small quark 
masses. We found that, for enough large spatial 
size 3.0fm), the mass splitting of the ground 
and first excited states of the nucleon is indepen- 
dent of the quark mass. This is consistent with 
the experimental mass splittings of baryons and 
mesons for wide range of quark masses By 
removing the finite volume effect, the mass of the 
Roper resonance becomes considerably reduced 
and the level switching between the Roper and 
the negative parity resonance may take place for 
sufficiently small quark masses. 

Needless to say, further study is required to 
check other systematic error stemming from a dis- 
cretization artifact. Additive simulations at dif- 
ferent lattice spacings are under way. 

We appreciate A. Nakamura and C. Nonaka for 
helping us develop codes for our simulations from 
their open source codes (Lattice Tool Kit 10 ). 
This work is supported by the Supercomputer 



Project No. 85 (FY2002) and No. 102 (FY2003) 
of High Energy Accelerator Research Organiza- 
tion (KEK). S.S. thanks for the supper by JSPS 
Grand- in- Aid for Encouragement of Young Scien- 
tists (No.15740137). 

REFERENCES 

1. For a recent review, S. Capstick and W. 
Roberts, nucl-th/0008028. 

2. S. Sasaki, T. Blum and S. Ohta, Phys. Rev. 
D65 (2002) 074503; S. Sasaki Nucl. Phys. 
Proc. Suppl. 83 (2000) 206; hep-ph/0004252. 

3. S. Sasaki, Prog. Theor. Phys. Suppl. 151 
(2003); nucl-th/0305014. 

4. S. Sasaki, K. Sasaki, T. Hatsuda and 
M. Asakawa, hep-lat/0209059. 

5. Y. Nakahara, M. Asakawa and T. Hatsuda, 
Phys. Rev. D60 (1999) 091503; M. Asakawa, 
T. Hatsuda and Y. Nakahara, Prog. Part. 
Nucl. Phys. 46 (2001) 459. 

6. T. Yamazaki et al. [CP-PACS Collaboration], 
Phys. Rev. D65 (2002) 014501. 

7. M. Fukugita et al, Phys. Lett. B294 (1992) 
380. 

8. M. Gockeler et al. [QCDSF-UKQCD-LHPC 
Collaboration], Phys. Lett. B532 (2002) 63. 

9. S. J. Dong et al, hep-ph/0306199. 

10. S. Choe et al, Nucl. Phys. Proc. Suppl. 106 
(2002) 1037. 



